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^ ' Abstract. In this paper, we consider operator realizations of quadratic algebras generated 

by second-order superintegrable systems in 2D. At least one such realization is given for each 
set of Stackcl equivalent systems for both degenerate and nondegenerate systems. In almost 
1/^ , all cases, the models can be used to determine the quantization of energy and eigenvalues 

for integrals associated with separation of variables in the original system. 
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^: 

1 Introduction 

^ , In his 1968 book Lie theory and special functions, W. Miller Jr. used function space realizations 

\ of Lie algebras to establish a fundamental relationship between Lie groups and certain special 

■ functions including Bessel functions, hyper geometric functions and confluent hypergeometric 

. functions [1]. It was further shown how the algebra relations can be used to identify special 

Tjj- \ function identities. 

' In this paper, we will apply these methods to the study of the representation theory for 

quadratic algebras generated by second-order superintegrable systems in 2D and their associated 
special functions. We would like to consider irreducible function space representations and so we 
restrict the Hamiltonian to a constant times the identity and construct difference or differential 
^ ■ operator realizations for the elements of the algebra. We call such operator realizations models. 

, A classical or quantum Hamiltonian on an m dimensional Riemannian manifold, with met- 

ric g^^ , given respectively by 

-H = \g'^PjPk + V{x^,X2), H = -^d^^^g^'^d,, + V{x^,X2) 

is called superintegrable if it admits 2m — 1 integrals of motion. For classical systems, we require 
that the integrals be functionally independent and, for quantum systems, we require that that 
they be algebraically independent within a Jordan algebra generated by Xj , dx^ and the identity. 
If both of the integrals are polynomial in the momenta in the classical case and as differential 
operators in the quantum case, we call the system nth-order superintegrable, where n is the 
maximal order of a minimal generating set of integrals. 

The study of superintegrability was pioneered by Smorodinsky, Winternitz and collaborators 
in the study of multiseparable systems on real Euclidean space [21 [3]. They identified all four 
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mulitseparable potentials on real Euclidean space, since named the Smorodinsky-Winternitz po- 
tentials. Later, W Miller Jr. with collaborators including E. Kalnins, J. Kress, and G. Pogosyan 
published a series of papers which classified all second-order superintegrable systems in 2D (see 
e.g. m m [6] and references therein). In these paper, the authors proved that classical and 
quantum second-order superintegrable systems are in one to one correspondence and that the 
potentials are invariant with respect to scaling so that it is possible to normalize the — /i^/2 term 
to be 1. Thus, for the remainder of the paper we will only consider quantum systems, except to 
make slight observations where the classical systems differ, and chose this normalization. 

It was further shown that any second-order superintegrable system in 2D can be related, 
via the Stackel transform, to a superintegrable system on a space of constant curvature and 
a complete list was given of all second-order superintegrable systems on 2D Euclidean space, 
£^2,0) and on the two sphere, S2,c- Since we will not use the Stackel transform explicitly in this 
paper, we refer the reader to (3 [8] and references therein for a complete exposition. We only 
note here that the Stackel transform is a mapping between Hamiltonian systems, possibly on 
different manifolds, which preserves superintegrability and the algebra structure of the integrals 
up to a permutation of the parameters and the energy. Thus, we can classify superintegrable 
systems based on the structure of their symmetry algebra. Such classifications have been worked 
out directly in [9l [10] and via the Stackel transform [12] . 

Most importantly for this paper, it was proved that the algebra generated by the constants of 
the motion for a second-order superintegrable system closes to form a quadratic algebra. That is, 
suppose H is second-order superintegrable with second-order integrals of the motion Li and L2. 
The integrals Li and L2 will not commute and we denote their commutator by R. The algebra 
to be considered is then the associative algebra generated by 

A = {Li,L2,H,R=[Li,L2]}. 

Such an algebra is called a quadratic algebra if the commutator of any two elements can be 
written as at most a quadratic polynomial in the generators. Further, since the four generators 
can not be independent in the Jordan algebra generated by the operators xj and dxj and the 
identity, there will be a polynomial relation between them. 
The structure relations for our algebra are hence given by 

[Li,R]=P{Li,L2,H), [L2,R] = Q{Li,L2,H), R" = S{Li, L2, H), (1.1) 

where P and Q are at most quadratic polynomials and S is at most cubic. We note that in 
the classical case, we will instead have a quadratic Poisson algebra and the highest order terms 
of the structure relations will be the same as in to the quantum case. If the system admits 
a first-order integral X, we call the system degenerate. In this case, the potential depends only 
on a single parameter and R can be expressed in terms of the basis X, Li, L2 and H [T3] . 

The study of the the algebras generated by superintegrable systems and in particular their 
representation theory has been a subject of recent study (see e.g. [HI [151 [TH [T71 [S] ) . The main 
advantage of the method of function space realization is that it can be used to find the eigenvalues 
of operators other than the Hamiltonian and to compute inter-basis expansion coefficients for 
the wave functions of the Hamiltonian. Further, as in the case of Lie algebras, the representation 
theory for quadratic algebras, and polynomial algebras more generally, seems to be intimately 
connected with special functions and their identities. 

This paper is divided up into four sections. In Section [21 we give a model which realizes 
the quadratic algebra associated to the singular isotropic oscillator. El or Smorodinsky-Winter- 
nitz i. A function space representation is given including the normalization and the weight 
function as well as the eigenfunctions and eigenvalues for each of the operators associated with 
separation of variables in the original system. We describe how the energy values for the system 
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and eigenvalues of the operators associated to separation of variables are quantized in a finite 
dimensional representation. In Sections [3] and HI we give at least one model for the algebras 
associated with a representative of each Stackel equivalence class of nondegenerate and degen- 
erate systems respectively. In these sections, the model is given and we identify possible finite 
dimensional representations including those are associated with quantized values of the energy. 
In Section [5] there is a brief conclusion with possible future developments. 

2 Exposition of a model: the singular isotropic oscillator 

In this section, we shall consider the quadratic algebra associated with the system El/Smoro- 
dinsky-Winternitz i and construct an irreducible function space representation for algebra. The 
Hamiltonian for this system, the singular isotropic oscillator, is 



y2 



The remaining generators of the symmetry algebra are. 

The algebra relations are, recall R = [Li, L2] , 

[R, Li] = 8Ll - 8HL1 + 16uj'^L2 - 8lo^, (2.1) 
[R, L2] = 8HL2 - 8{Li, L2} + (12 - 16a'^)H+ (iGa^ + 166^ - 24) Li, (2.2) 

R^ = 8H{Li,L2} - -{Li,Li,L2} + 16uj^Ll + 16{a^ - l)H^ + ll6a^ + 166^ - — \ if 

_ /'32a2 - HLi - ^^L2 - — {48aH'' - 48a^ - 486^ + 29) . (2.3) 

Y 3 / 3 3 

In order for the representation to be irreducible, the Hamiltonian will be restricted to a con- 
stant and so, in the model, the action of the Hamiltonian is represented by a constant E, 
i.e. H = E. Now, suppose that Li is diagonalized by monomials, that is, it is of the form 
Li = litdt + I0 where h and Iq are some constants to be determined later. Suppose also that the 
action of L2 on the basis of monomials is in the form of a three-term recurrence relation. Prom 
the algebra relations, it is straightforward to see that L2 cannot be modeled by a first-order 
operator or else the relations would be a Lie algebra. Thus, we take the Ansatz of a second-order 
differential operator with polynomial coefficients. 

The algebra relations (|2.ip - (j2.3p are realized by the following operators, H = E and 



Li = -Autdt + E- 2u){l + b), 



2 V cj 2 



2t((E-4uj-2Lob)^-4oo'^a'^) 2E(l + b) ^,0 
H ^ + 46^ + 86 + 5. 

As was proposed in the Ansatz, Li is diagonalized by monomials in t and further the opera- 
tors Li and L2 act as automorphisms on polynomials in t. Thus, we take the set {t^ | A: = 0, . . .} 
to be the basis for our representation. The three-term recurrence relation generated by the 
action of L2 on monomials is 
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where 



2{E - 4a; - 2aw - 2bu - 4na;)(S - 4u; + 2auj - 2huj - Anu) 

Cn 



-IGn^w + {AE - 16a;(l + h))n - 2E{1 + b) + uj{4b^ + 86 + 5) 



Cn— l,n 



2u; 

n(?i + b) 



Note that for arbitrary parameters, C_i^o vanishes whereas Cn+i,n in general does not vanish 
for any value of n and so the so the representation is infinite dimensional, bounded below. If we 
make the assumption that Cn+i,n is zero for some finite integer, then we obtain quantization 
conditions for the parameters and a finite representation of dimension m. In particular, we 
chose to solve Cm,m-i = for the energy value E to obtain 

E = -2oj{2m + a + b), m G N. (2.5) 

The potential and hence the algebra relations are symmetric under the transformation b — t- —b 
and a —7- —a and so we could have taken the opposite sign for either a or 6 for the construction 
of the model and hence in the resulting quantization of the energy values. 

The existence of a three-term recursion formula as in (|2.4p indicates the existence of raising 
and lowering operators. Indeed, they are given by the following relations 

^"^^ + £-|^^+2^^^-^=*^*+(^ + ^)^^' 

-r R Li E ^ 1 

= Lo k H T^Li 

= 64*^52 _ 3'^{E-2ub-6u) ^^^^ ^ 4{{E - 2ujb - Au)^ - Aa^u) 

Here we note that if L\ and L2 are self-adjoint, A and will only be mutual adjoints if uj 
is real, which is reflective of the physical fact that the potential will be attracting in that case. 
The commutation relations of the raising and lowering operators can either be determined from 
the quadratic algebra (|2.1P " (I2.3P or directly from the model. They are 

[Li,y4] = AijjA, \Lx,A^\ = -AijjA\ 

[A, At] = + - A [E^ _ 4a;2(,2 + ft^ + 2)) - ^1^1^, 

OJ-^ UJ-^ (jj-^ UJ 

We return to finite dimensional model where the basis of eigenvectors for Li, is given by 
{i;^n(i) = I n = 0, . . . , m — 1}. Let us assume the existence of an inner product for which L\ 
and L2 are self-adjoint, or practically, we assume A and A^ are mutual adjoints. As mentioned 
above, this is equivalent to L\ and L2 being self-adjoint and the constants a, 6, and uj being real. 
With this inner product, we can find the normalization for our eigenvectors using the raising 
and lowering operators. That is, we assume 



n-l\ 



and = 1 to obtain the recursion relation as 



2 64(m — n)(m — n + a) 2 

= — - K, 



n{b + n) 



n— 1 J 
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so that 



16'' 



' (-m)„(-m - g), 
n!(6)n 



which is real so long as a,b > 0. 

From these normalization coefficients, we can find a reproducing kernel for this Hilbert 
space which lies in the Hilbert space. It has the characteristic {6{ts), f{t)) = /(s) and is given 
by 4'n{t)4>n{s) = knits)^ which is exactly the hypergeometric polynomial 



5{t,s) 



2F1 



-m, —m — a 
b 



ts 



Finally, it is possible to construct an explicit function space representation for the inner 
product which will make Li and L2 self-adjoint. We assume the inner product is of the form 
{f{t),g{t)) = f{t)g{t)p{t,t)dtdi where 7 is a path to be determined later. We can determine 
the weight function using the following relation, 



(Li/,g)= / {Lif{t))g{t)p{t,t)dtdt= / fit){Lig{t))p{t,t)dtdt = {f,Lig). 
J ^ J ^ 

Similarly, we require {Af,g) = {f,A^g). These two integral equations can be transformed by 
integration by parts into differential equation for the weight function whose solutions are 



p{tt) = ci 2F1 



l + m, m+l + o 
1 - b 



tt +C2{ttf 2F1 



1 + 771,-1-6, m + 1 + a + b 
1 + 6 



tt 



Finally, we note that we can also diagonalize the operator L2 in the model. The solutions of 
the eigenvalue equation for L2, {L2 — X)ipx = 0, are hypergeometric functions, 2-^1 's [E]. If we 
assume the model is finite dimension and restrict E to the values in (|'2.5|) . then L2 becomes 



t{8t - 1) 



-0? 



{8t + 1)((2777 + a - 3)8t -6-1) 



dt 



+ 32(m - l)(r77 - 1 + a)t - 4m{b + 1) - 26a - 2a + 26 + 



For a finite dimensional irreducible representation, the hypergeometric series must terminate and 
we obtain a quantization relation on the eigenvalues A = — 3/2— 26— 2a— 4/c— 26a— 46A:— 4afc— 4/i;^ 
and the basis functions become, for A: = l,...,777 — 1, 



Mt) = k{8t + 1) 



m—l—k 



2F1 



—k, —a — k 
1 + 6 

lb, a 



8t 



k{8t + l 



^m-lpb,a / 1 — 8t 
^' ^ l + 8t 



Here the l^s are normalization constants and P^' is a Jacobi polynomial. 

This model is an interesting example of the simplicity of differential models. We can di- 
rectly compute the eigenvalues for all the operators and the normalizations of the basis vectors. 
Also, the kernel function and integral representation of the inner product give useful tools to 
find expansion formulas and possible generating functions for the original quantum mechanical 
problem. 

In the remainder of the paper, we exhibit a function space realization for a representative of 
each of the Stackel equivalence classes of second-order superintegrable systems in 2D. 
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3 Models of quadratic algebras for non-degnerate second-order 
superintegrable systems 

We will make use of the following conventions. On Euclidean space, the generators of the Killing 
vectors are, 

Px = dx, Py = dy, M = xdy - ydx- 

We also define the operator p± = dx ± idy. 

The Laplacian on -E'2,c in Cartesian coordinates is 

In complex coordinates, z = x + iy, z = x — iy, we have 

A = d,d^. 

Or, we can take a real form in lightlike coordinates, v = x, C = iy, to obtain the wave operator. 

For systems on the two sphere, we use the coordinates of the standard embedding of the 
sphere into 3 dimensional Euclidean space. We denote these si, S2, S3 such that sf + S2 + s'^ = 1. 
The basis for the Killing vectors is 

J. = 6^^- Vs.- 

The Laplacian on 52,c is 
1=1 

The quantum algebra structures often have symmetrized terms. We define these by {a,b} = 
ab+ba and {a, 6, c} = abc+acb+bac+bca + cab + cba. Also used is the permutation sign function 
€ijk and eijki, the completely skew-symmetric tensor on three or four variables respectively. 

In the following sections, we exhibit irreducible representations of the quadratic algebras, 
defined by relations (jl.ip . using both differential and difference operator realizations of the 
operators acting on function spaces. Because the Hamiltonian operator commutes with all 
of operators, it must be a constant for any irreducible representation and so we restrict the 
Hamiltonian to a constant energy H = E. We focus on models which will diagonalize operators 
associated with separation of variables, since these are of immediate interest for explicit solution 
of the physical system though additional models are given as well. Most of the models were 
obtained by assuming a differential operator Ansatz except in the case of S9 which is realized 
as a difference operator. In this case, the operators were derived from the abstract structural 
relations though it is interesting to note that they also could have been obtained through the 
quantization of a model for the Poisson algebra for the associated classical system [20^ I21j . 
Unless otherwise noted, these models were first exhibited in |22| . 

We classify the nondegenerate systems in 2D by their Stackel equivalence classes and identify 
them by the leading order terms of the functional relation, as seen in the accompanying Table [TJ 
This classification comes from [11] and the nomenclature of the Ei's comes from [3]. In this 
section, we describe the symmetry operators, quadratic algebras and at least one model for each 
of the equivalence classes. 
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Table 1. Stackel equivalence classes of non-degenerate systems in 2 dimensions. 



Leading terms of Casimir relation 


System 


Operator models 


Ll +f{auH)Ll 
L\ +f{ai,H)L^L2 
Ll +0 
LIL2 +f{ai,H)Ll 
Lf L2 + 
LiL2iLi+L2) + f{ai,H)L,L2 
+f{ai,H)LiL2 


E2, SI 
E9, ElO 
E15 

El, E16, S2, S4 
E7, E8, E17, E19 
S7, SB, S9 
E3, Ell, E20 


Differential 

Differential 

Differential 

Differential 

Differential 

Difference 

Differential 



3.1 E2: Smorodinsky— Winternitz ii 

The Hamiltonian is on real Euclidean space for real constants u, b, c 
H = A- u?{Ax^ + y^) + bx + 



y 

A basis for the symmetry operators is given by H and 

Ll = pI - Auj'^x'^ + bx, L2 = ^{M,py} - Q - x^^^ + (i " ''^) J- 

The symmetry algebra relations are, recall R = [^1,^2], 

[Li,R] = -2bH + 16Lj'^L2 + 2bLi, (3.1) 

[L2,R] = 8L1H - 6Ll - 2H^ + 2bL2 - Suj"^ {I - c^) , (3.2) 
R^ = 4L? + ALiH"^ - SLjH + IGw^L^ - 'ibL2H + 26{Li, L2} 

+ 16a;2(3 - c^)Li - 32uj'^H - 6^(1 - c^). (3.3) 

The algebra relations ()3.ip - ()3.3p are realized hy H = E and the following operators, 



1^ 

16^2 



L2 = 32ujt^d^ + ( ( 16{E - 6w) + -^"j i2 _ At _ i ) dt 



't 

,2 



^ {l&Euj'^ - 64^3 _ ^2)2 _ (32a;3g)2^ ^ IQEuj"^ - 2,2u'^ - b' 



128tj5 I28tj4 

Here, the eigenfunctions of Li are monomials and the action of L2 on this basis can be represented 
as a three-term recurrence relation given by 

where 



n+l,n 



16a;2 



_ n r - K16L;w2 - 32^^ _ 52 _ 64a;2„) 



128^^4 

Notice, that C_i^o vanishes and so our representation is bounded below. Further, if there 
exists some n so that Cn+i,n also vanishes then the representation becomes finite dimensional. 
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Conversely, we can assume that the representation is finite dimensional to obtain quantization 
conditions on the energy. If we assume that representation space is m dimensional and spanned 
by the monomials {i" \ n = 0, . . . ,m — 1}, then the restriction Cm,m-i = gives the quantization 
of the energy values 

E = Aujim + 2ec) H ^, m G N. 

Here, we can take either e = ±1 in the energy which is consistent with the Hamiltonian depending 
only (?. 

Finally, we note that we can define raising and lowering operators in this model as, 
R hLi bE _ 1 

. , R bLi bE 

= t {Uujt^d\ ^ ^^-5 -tdt + [X^Eu? - UiJ -\?) - (32a;3c) j . 

The raising and lowering operators obey the following commutation relations 
[Li , ^] = -4wA, [Li ,A^\= Au)A^ , 

3.2 ElO 

The Hamiltonian for this system is given by, with z = x + iy, z = x — iy 

3_2\ , 1-3 



i? = A + az + /3(^z- -z^'j + 7 l^zz - -z 

A basis for its symmetry operators is given by H and 
Li=pi+ iz^ + 2/3z, 

3 

L2 = 2z{M, } + p+ - 4/32;z - 7zz^ - 2/3z^ - -72"^ + 72^ + az^ + laz. 
The algebra relations are given by 

[i?, Li] = -327L1 - 32/32, ^ gg^2 _ i28aLi + 327L2 + 64/3i? + 32a^ 

R^ = 64L? + 327{Li, L2} - 128aLf - 647^^^ _ US^HLi + 64/3^^2 + 6Aa'^Li 
- USPaH - 2567^. 

This algebra can be transformed into the Lie algebra s/2 by using the following invertible trans- 
formation 

In this basis, the algebra relations reduce to a Lie algebra 

= -327A'i, IJi.ffj] =327/1-2, [Ki,ff2l = R, (3.4) 
= 327{if,. Ay - 647i/^ - 1M21±«^ _ 64(/i° + 4^* + „ W + 2./i^0 ^ 
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The algebra relations (|3.4p . ()3.5p are realized hy H = E and the following operators, 

Ki = lQ-idt, K2=t dt+ \ 

~(7^g + a/37 + p"^] 

The operator R is diagonalized by monomials t" and in this basis Ki and K2 act as lowering 
and raising operators respectively. Note that Ki annihilates constants and so the representation 
is bounded below. On the other hand, for a finite dimensional representation, we require that 
there exist some integer m such that 



i?=[Li,L2] = 167 2t9t+ 1 + 



i^2t--^=lm-l+ll+ - - — ' - ) )t- = 0. 



-7(7^^ + a/37 + 1^^ 
273 



This leads to the quantization condition on the energy 

a/3 /33 

E = 2iJ^m meN. (3.6) 

7 7 

Notice that in order to obtain a real energy value we require that all the parameters be real and 
that 7 < 0. 

It is also possible to diagonalize a linear combination of Ki and K2 - For example, the solutions 
of the eigenfunction equation 

{Ki + - A)^ = 



are 



\I/ = (167 + t^) ir^ exp ( — ^ arctan 



4^7 V4V7 

which, for finite dimensional representations where E is restricted to the value in (13. 6|) . give 
quantization conditions on the eigenvalues A. That is, if we require that \I' be a polynomial in t 
of degree less that m we obtain a complete set of eigenfunctions given by 

{Ki+ K2)-^n = Xn'^n, ?1 = 0, . . . , 771 - 1 , 

^'n= (4^/^ + t)"(4V^-t)'"-"-\ A„ = 4V^(?n-2n-l). 

Finally, we can return to the original basis of Li and L2. In the model, 

/32 

Li = IQ-idt , 

7 

L2 = 256^8^ + (^t^ + 32a + dt 

( ^(7^g + a/37 + /33) ^^ _ 2/37^g + a^7^ + 4a/327 + 3/3^ 

+ V 273 J* 73 

3.3 E15 

Here 



H = A + h{z), 
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where the potential is an arbitrary function of z. A basis for the symmetry operators is 

/dh 

The only nonzero algebra relation is [Li,L2] = iLi. This system is unique among all 2 dimen- 
sional superintegrable systems in that the symmetry operators are not functionally linearly in- 
dependent and do not correspond to multiseparability. The only separable system is determined 
by diagonalizing Li, essentially 2, and this coordinate system is not orthogonal. A model is 

d r -.d . , 

Li = — + a, L2 = ^^— -|- lat, 
at at 

but the irreducible representations of the algebra yield no spectral information about H. Again, 
since this algebra is a Lie algebra, the above model is derivative of one found in [l] . 

3.4 El, Smorodinsky— Winternitz i 

This is the system considered in Section [2] and for completeness we recall the results here. The 
Hamiltonian for the system is 



4 



The remaining generators of the symmetry algebra are 



,2 , 1/4 - a' , 02 , _ ,,2 , (1/4 - a^)y' , (1/4 - b^)x^ 



The algebra relations are 

[R, Li] = 8Ll - 8HL1 + 16uj'^L2 - Suj"^, (3.7) 

[R,L2] = 8HL2 - 8{Li,L2} + (12 - 16a^)H+ (iGa^ + 166^ - 24) Li, (3.8) 

i?2 = 8H{Li,L2} - ^{Li,Li,L2} + IGw^L^ + 16(0^ -l)H^+ (^16a^ + 166^ - Lj 

( n 200\ 176^^2 4a;2 , ox . . 

- 32a^ HLx L2 {ma%^ - ^a^ - m? ^ 29) . (3.9) 

V 3 / 3 3 



The algebra relations (|3.7p - (j3.9p are realized hy H = E, and the following operators 
Li = -Aujtdt + E- 2uj{l + b), 

L, = i*(8t + If a! - + _ 111] a, 

2t((E -2ub)^ -Auj^a^) 2E(l + b) ,0 
H TT-^ ^ + 46^ + 8b + 5. 

LO'' UJ 

As described in the previous section, the eigenvalues for Li are monomial and the eigen- 
functions for L2 are hypergeometric functions which reduce to Jacobi polynomials for finite 
representations. The representation become finite dimensional under the quantization of energy 
E = -2u{2m + a + b), (f23D . 
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3.5 E8 

The Hamiltonian is 

H = A + ^ + + jzz. 

Z"^ z 

A basis for the symmetry operators is given by H and 



z z z 

The algebra relations are 
2 



Li] = 8Li^ + 32a7, [i?, L2] = -8{Li, L2} + 86i7 - 16(a + l)Li, (3.10) 

8 / 176\ 

= --{Li, Li, L2} - \ lQa + — \ lI + 16aH^ - 6Aa-fL2 + IQpLiH 

The algebra relations (jS.lOp . (jS.lip are realized hy H = E and the following operators 

= = -4(«^ - m + ((^ - 8) * + ^) a. - (1 + «. 

Here, the eigenfunctions of L2 are hyper geometric functions which restrict to Jacobi polynomials 
under quantization of eigenvalues 

L2^ =A^ A = -4n2 + ('^^ - ^^"^ ^ _ + 4a + 4/3^^ - /g^ 



The eigenfunctions are given by 



4^/^ 4^' 4^/^ 4^' 



where is a normalization constant. The action of Li in this model is via multiplication by 
the variable t and gives a three-term recurrence formula. Note that for some quantization of the 
energy the eigenfunctions become singular for all n > m. This occurs under the quantization 
condition m E N with energy eigenvalues 

E = 2\f^ I 2m + 2 ± , m G N. 

3.6 S9: the generic system in 2D 

The Hamiltonian is 



1 a2 1 _ 52 12 
H = As2 + + + 4 

"1 ^2 



A basis for the symmetry operators is 



i2 = J/ + l.-oMi|+fT-''^)4 



3 

H = Li + L2 + L3 + --a^ 
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The structure equations can be put in the symmetric form using the following identifications 

1 2 -^2 ^ i2 

ai = --C , a2 = - -a , 03 = --b , 

[L„R] = eijk (4{L„ Lk} - 4{Li, Lj} - (8 + 16aj)Lj + (8 + 16afc)Lfc + 8(0^ - 0^)) , (3.12) 

i?2 = ^{Li, L2, L3} - (16ai + 12)L? - (16a2 + 12)L^ - (IGas + 12)L^ 

52 11 
+ — ({Li, L2} + {L2, L3} + {L3, Li}) + 3 (16 + 176ai)Li + -(16 + 176a2)L2 

1 32 
+ -(16 + 17603)^3 + — (ai + 02 + 03) + 48(aia2 + 0203 + asoi) + 64aia2a3- (3.13) 

We can obtain Li in the model by using the following difference operator, based upon the 
Wilson polynomial algebra. This model is unique in that the energy values E and eigenvalues 
of the operators were used in determining the model. We define the coefficients a, /3, 7 and 6 
as 

a+c+1 ^ a+c+1 a-c+1 ^ a+c-l 
" = ^ P = ^ , 7 = ^ , d = + 6 + ^ + 2, 

and use difference operators 

T^F{t) = F{t + A), r = -(ri/2_r-V2), 

T* = ^ [(a + t){(3 + t)(7 + t)(5 + t)ri/2 _ (a _ t)(;3 - t){^ - t){6 - t)T-V2] . 

The algebra relations (|3.12p . (|3.13p are realized by the following operators 

L3 = -4*2 + + c^, Li = -4r*T - 2(a + 1)(6 + 1) + ^, 

= E = -^(4/i + 2a + 26 + 2c + 5)(4/i + 2a + 26 + 2c + 3) + ^ - - 6^ - c^. 

The model realizes the algebra relations for arbitrary complex fi and restricts to a finite 
dimensional irreducible representation when /i = m G N. In this model, we obtain spectral 
resolution of L3 with delta functions as eigenfunctions. The eigenfunctions of Li are Racah 
polynomials in the finite dimensional case and Wilson polynomials for the infinite dimensional, 
bounded below case. This model was first published in [23] where the model was worked out in 
full generality including the normalizations and the weight functions. It has also been recently 
been extended to the 3D analog in |16j. 



3.7 E20: Smorodinksy— Winternitz iii 

The Hamiltonian is 



\/ x2 + ■u2 _|_ 2; \/ J x2 + %fl — X 

H = A+ ' I a + ^ + 7 ^ ^ ^ - 

Vx2+y2 1 ^a;2_|_y2 ^3;2_|_y2 



In parabolic coordinates (x', y'), with x' = ysP^+y^ + x, y' = \J x^ + ip' — x, the Hamiltonian 
can be written as 
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The integrals are 

^2 = ^72^ {-^'y'{d!' + dy') + + y")d.'dy, - 4ax'y' + 2{x'' - y") (y'/3 - x'7)) • 
The algebra relations are given by, 

[R,Li] = -4L2H + 16/3^, [R,L2]=iLiH -8{P^ -j^), (3.14) 
= ALIh + ALIh + AH^ - IGa'^H + 16(7^ - /3^)Li - 32/37L2 - 32a'^{/3^ + j^). (3.15) 

Notice, these restrict to a Lie algebra on a constant energy surface, H = E, and so the model 
described below is not new but, for example, can be derived from those given in [1]. 

The algebra relations ()3.14p . ()3.15p are realized hy H = E and the following operators. 



q2 

Li = -2\fEtdt - \/E + 2a + A--, 

E 



= - \Et-d, + 2d, - \tE + V^at + + ^) t - 4 



/37 



E 



Here L\ is diagonalized by monomials and the action of L2 on monomials is given by the 
following three-term recursion formula 

with 

a + 7^ 4/37 

Cn+l,n = —{n + l)E -\ 7= — , Cn,n = —^, Cn-l,n = 2n. 

\/E \/E -c/ 



Notice, that C-ifi vanishes for arbitrary parameters and so our representation is bounded 
below. On the other hand, if we require that the representation be finite dimensional, say of 
dimension m, we obtain the restriction 

r- 0^ 7^ 

-mE + 2VEa + -= + ^ = 0, m G N, (3.16) 
VE ^/E 

which gives the quantized restrictions on the energy. 
The eigenfunctions for L2 are given by 

^ = {Et^ - |^^^_i^arctanh (^^^ ^ , 

where L2'^ = A^. If we assume that the model is finite dimensional, i.e. that the energy value E 
satisfies (|3.16p . then the eigenvalues are quantized as 

A = A„ = ^ + (m-l - 2n)/E 
E 

and there exists a complete set of eigenfunctions satisfying 
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This model admits raising and lowering operators, 
^ = L, + ^-^ = 45„ 

2VE E 
which satisfy the commutation relations 

[A, Li] = -2VeA, [A\Li] = 2VeA\ 



4 2D degenerate systems 

Next, we consider degenerate systems whose symmetry algebra includes a first-order integral, X. 
As described above, this requires that the potential depend on only one parameter, not including 
the trivial additive constant. In these systems, the symmetry algebra is defined by the operators 
H, Li, L2 and X. It is always possible to rewrite the commutator R = [Li,L2] as a polynomial 
in the other operators of maximal degree 3 in X and 1 in Li, XLi, L2, XL2, H and XH. In 
these systems, the commutation relations are in terms of X instead of R. That is, the defining 
relations are 

[Li,X]=Pi{Li,Li,X^,X,H), [Li,X]=P2{Li,Li,X^,X,H), 
[Li,L2] = Q[X^,XLi,XL2,XH,Li,L2,H,X). 

Here, the Pj's and Q are linear in the arguments. Furthermore, the functional relation is no 
longer in terms of R^ but instead a fourth-order identity. 

We begin with the table of the equivalence classes of degenerate systems; there are exactly 6 
degenerate systems in 2 dimensions. 



Table 2. Stackel equivalence classes of degenerate systems which admit a Killing vector X, in 2D. 



Leading order terms 


System 


Operator models 


+ L1L2 + AX'^ 


E3, E18 


Differential 


X^ + L1L2 


S3, S6 


Differential and Difference 


X^ + X^Li +Ll +0 


E12, E14 


Differential 


+ X'^Li + + ALi 


E6, S5 


Differential 


X^ +Ll 


E5 


Differential 


+X'^Li +L2 


E4, E13 


Differential 



4.1 E18: the Coulomb system in 2D 

This system is defined by the Hamiltonian 



^3,2 _|_ y2 

A basis for the symmetry operators is formed hy H, X = M and 
Li = 7t{M,p4 - ^ r^^^ , L2 = -{M,py} 



2' '"^ 2^x2 + 2/2' 2^ '"^^ 2^2:2 + y2 
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The symmetry algebra is defined by the following relations 

[Li,X]=L2, [L2,X] = -Li, [LuL2] = HX, (4.1) 

Ll + Ll-HX^ + ^^^-^ = Q. (4.2) 

We can change basis so that the algebra is in the standard form of the Lie algebra 5/2- With 
the substitutions A = Li + iL2, = Li — iL2 the algebra relations ()4.ip . (j4.2p become, 

[A,X] = -iA, [A\X] = iA\ [A,A^ = 2iHX, (4.3) 

{A,A^}-2HX^ + ^-^ = 0. (4.4) 
The Lie algebra defined by (j4.3p . (j4.4p is realized hy H = E and the following operators 

X = -i(tdt + l + ^] , A = ^dt, A^ = -4t(tdt + i+ " 



2 2VEj' 4 " \ VE 

The operator X is diagonalized by monomials and for arbitrary E and a the representa- 
tion is infinite and bounded below. Finite dimensional representation occur if A^^ annihilates 
a monomial. That is 

^t^m-i ^ 4(^m - aVE) ^ ^ 
E 

when the energy value E takes the values 

E = m G N. 

Again, since this algebra is a Lie algebra the model employed is not new. 

4.2 S3 

This system is defined by the Hamiltonian 

4 

A basis for the symmetry operators is formed hy H, X = J3 and 

Li = Jf + — 2"^-^) L2 = -{J1J2 + J2J1) — — 2^ • 

S3 I S3 

The algebra relations are given by 

[Li,X] = 2L2, [L2,X] = -X^ -2Li + H + a^ (4.5) 
[Li, L2] = -{Li,X} + (2a2 - 1)X, (4.6) 

^{Li, X, X} - HL^ + Ll + Ll- {a" " ^) " + ^) ^1 

The algebra relations (j4.5p - (j4.7p are realized hy H = E and the following operators 
X = 2itdt + ico , 
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Li = t{t + lydt + {t + l)(cit + ci + 2co - 1)^4 

3cl E 19\ a'^ + d + E 1 



+ ^c? + ^-2coci + 3co-3ci-- + - + -jt+ ^ 

L2 = - t)dt - i {cit^ + 2 + 2co - ci) St 

/ 9 3c^ a2 S 19\ 

- W c? + ^ - 2coCi + 3co - 3ci - y + - + — j t, 

where, for compactness of the equations, we have chosen constants cq, ci as 



Co = a - IH , 

3 1 / 

ci = - (2a + \/4£; - 1) + |\/4a2 + 16a - UE - 13 - 4(a - 2)V4S - 1. 

In the model, X is diagonalized with monomials and both Li and L2 have a three-term recurrence 
relation in this basis given by 

= ^C„+i,„t"+i - iCn.,,nt''-\ (4.9) 

with 

^ 2 / N 2 3c^ a^ + 19 
<-^n+l,n = n + [ci - l)n + Ci + — 2coCi - 3ci + 3co ^ h 

-|- -£/ + q2 

Cn,n = 2n^ + 2nco + — -, C„_i,„ = n{n - ci + 2co + 1). 

Note that, for arbitrary parameters, C_i^o vanishes and so the model is bounded below. On the 
other hand, if we require that the representation space be finite dimensional, i.e. Cm,m-i = 0, 
we obtain the quantization condition 

E = -{m-af + ^, mEN. (4.10) 

As can be seen directly form the three-term recurrence formulas (|4.8p and (j4.9p . there exist 
raising and lowering operators given by 

= Li + iL2 + ^ (^X^ - E + ^ - a^^ , A = Li - iL2 + ^ (^X^ - E + ^ - 

The symmetry algebra of the raising and lowering operators is given by 
[A,X] = 2iA, [A\X] = -2iA\ 

{A, A^} - ^X^ + X^ (^E -a"^ + + ^{AE + + 8a + 3) (E + 40^ _ 8a + 3) = 0. 

There is also a difference operator model associated with the spectral resolution of L\. The 
operators are formed from the operators defined as T^f{t) = f{t + A;). When restricted to 
finite dimensional representations, i.e. the energy E takes the values as in (|4.10p . the algebra 
relations (I4.5p ~( j47f|l are realized by 



_ ^ (l/2-a-t)(m + a-l/2-t) ^i _ (1/2 - a + t)(m + a - 1/2 + t)) 
^ 2t 2t ' 
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(1 - 2t)(l/2 - a - t)(m + a - 1/2 - t) , 
-L'l — ; 

^ (l + 2t)(l/2-a + t)(m + a- l/2 + t)) ^_i 
4t 

The eigenfunctions of Li are delta functions and the eigenfunctions of L2 are dual Hahn polyno- 
mials. This model was first published in [20] where the model was worked out in full generality 
including the normalizations and the weight functions. 

4.3 E14 

The Hamiltonian is, with z = x + iy, z = x — iy, 
H = A + ^. 

A basis for the symmetry operators is formed by if, X = p_ and 

2 iz z 

The symmetry algebra is 

[X,Li] = iX2, [X,L2] = 2iL2, [Li,L2] = i{X,L2} + ^X, (4.11) 

L?-^{L2,x2} + aif-^x2 = 0. (4.12) 

The algebra relations (|4.11|) . (|4.12|) are realized hj H = E and the following operators 

a: = Li = idt, L2 = -t^dj - 2tdt + aEt^ - 

Eigenfunctions for L2 are Bessel functions while eigenfunctions of Li are exponentials both with 
continuous spectrum. This model has no obvious finite dimensional restrictions. 

4.4 E6 

The Hamiltonian is 

i/ = A + ^^. 



A basis for the symmetry operators is formed hy H, X = py and 
The algebra relations are given by 



[Li,X]=H-X', [L2,X]=2Li, [Li,L2] = {X,L2} + {l-2a)X, (4.13) 
Lj - HL2 - 2{Li,X} + 1^ + Q - A^ = 0. (4.14) 

The algebra relations (j4.13p . (|4.14p are realized hj H = E and the following operators 
X = VE{t^dt + ia + l)t + l), 

Li = -VE{t^d^ + ((2a + 3)t + 2)tdt + {a + ift + a + l) , 
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L2 = -t^dl - 2((a + l)t + l)dt - a - ^• 
Note that the operators X, Li and L2 act on monomials via three-term recurrence relations 

Xe = y/E{n + a + 1)*"+^ + Ve^, 

Li^ = -\/E{n + a + ife^^ - /E(2n + a + 1)^, 

Lst" = -{n'^ + {n + 2a) + a + ^ - Inf-^ . 

Thus, the model can be realized as a bounded below representation using the set of mono- 
mial {t" \n E N}. The representation will become finite dimensional, with dimension m if the 
parameter a is restricted to quantized values a = —m. Notice, a constant is quantized instead 
of the energy which does not make sense in the physical system, since the constants should be 
given by the system. However, under a Stackel transform the constant and the energy will be 
interchanged so our model is giving quantization levels of the energy of a Stackel equivalent 
system; in this case, S5. The eigenfunctions of Li are gauge equivalent to Laguerre polynomials 
with eigenvalues ^/E{2n — m + 1). The eigenfunctions of L2 are gauge equivalent to generalized 
Laguerre polynomials with eigenvalues — k'^ + k — ^. 

4.5 E5 

The Hamiltonian is 
H = A + ax. 

A basis for the symmetry operators is formed hy H, X = py and 

Li = PxPy + ^ay, L2 = ^{M,py} - ^ay^ 

The algebra relations are given by 

[Li,X] = -^, [L2,X]=Li, [Li,L2] = 2X^ - HX, (4.15) 
X'^ - HX'^ + lI + aL2 = Q. (4.16) 

The algebra relations (|4.15p . (|4.16p are realized hy H = E and the following operators, 

X = t, L^ = -^dt, L2 = -]ad^ + -f{E-t^). 

2 4 a 

Here the only nontrivial eigenfunction is that of L2 which gives solutions of the triconfluent 
Heun equation. This model gives no obvious finite dimensional restrictions. 

4.6 E4 

The Hamiltonian is 

H = A + a{x + iy). 
A basis for the symmetry operators is formed hy H, X = p^ and 

Li=pl + ax, L2 = ^{M,p+} + ^-^{x + iyf, 
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The algebra relations are given by 



[Li,X] = -a, [L2,X]=iX', [Li,L2] = -iX'' -iHX + i{Li,X} , (4.17) 

X"^ -- {Li,X, X} + 2HX'^ + H'^ + 4mL2 = 0. (4.18) 
3 

The algebra relations (|4.17p . (|4.18p are realized hy H = E and the following operators 

... r E ^ 2o iaH'^ - ^aH + iE"^ 

X = lat, Li = idt + — , L2 = -at dt H . 

2 4a 

This model gives no obvious finite dimensional restrictions and the basis for the representation 
can be chosen as exponential functions. 

5 Conclusion 

In this paper, we have demonstrated at least one function space realization for each class of 
Stackel equivalence superintegrable systems. We have shown that these quadratic algebras are 
related to various families of special functions and orthogonal polynomials. 

There are several possible directions that this analysis can be extended and the model de- 
scribed above can be further applied and analyzed. One possible direction, is to use the fact 
that there exist appropriately chosen limit processes which takes superintegrable system S9 on 
the sphere to all other superintegrable systems in 2D (see e.g. |24| I25|). It would be interesting 
to see how these limits are manifested at the level of the algebras and particularly the models. 
Also of interest would be to consider the restrictions from nondegenerate to degenerate systems 
and to study the effect of such a limit on the models. 

These methods can also be applied to the construction and analysis of quadratic algebras 
generated by superintegrable systems in higher dimensions. In particular, the quadratic algebras 
for all real nondegenerate 3D superintegrable systems were shown to contain various copies of 
quadratic algebras from two dimensional systems [26] . As was the case in [161 [27] , the models for 
the 2D system can be used as a basis for the models of the 3D system. It would be interesting 
to see how the algebra decomposed with respect to the subalgebras and whether that has any 
effect on the models obtained and possible special function relations. 
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